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' Abstract 

We determined the exact location of the shearless KAM curve in the quadratic map and 
numerically investigate the breakup thresholds of those curves in the entire parameter space. 
,_H , The breakup diagram reveals many sharp singularities like fractals on the reconnection 

^ ■ thresholds of the twin-chains with rational rotation numbers. 

1 Introduction 

m ■ 
O ' 

£ — , | In this paper we focus on the complicated phenomena which occur in the hamiltonian systems 

beyond the twist condition. Because of the violation of the twist condition |l], the KAM the- 
orem and the Poincare-Birkhoff theorem do not hold near the shearless KAM curve where 
the rotation number takes an extremum value. When a perturbation is added, very complicated 
phase space structures are observed in the vicinity of the shearless KAM curve: twin resonant 
chains appear on both sides of the shearless KAM curve, and at a certain parameter value the 
separatrices of the twin-chains merge. These phenomena, so-called reconnections of twin-chains, 
' are generic ones in many nontwist systems f|]-f§. 

The quadratic map, which is the simplest example of a nontwist system, has been studied 
from the viewpoint of the reconnection of twin-chains, and the breakup mechanism of the shear- 
less KAM curve has been studied to a certain extent^, ||, ||. However the detailed phase space 
structures and the breakup condition of shearless KAM curves in the entire parameter space 
have not yet been elucidated completely. 

In this paper, carrying out with the quadratic map, we derive the exact location of the shear- 
less KAM curve and numerically determine the breakup threshold of shearless KAM curves in 
the entire parameter space. 

The quadratic map is the following two-dimensional area-preserving map (I n , 9 n ) \— * (I n +i, #n+i) 




- Ksm(6 n ), 
On + fn{I n +i), (mod 2tt) 



(1) 



/ At (/) = 27T^-/ 2 



where K represents the strength of the perturbation, and \jl is the rotation number of the 
shearless KAM curve in the integrable case, i.e.,K = 0. Because the variable 9 is 27r-periodic, 
fi is also periodic with period 1. Thus it is enough for us to consider the parameter region 
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-0.5 < n < 0.5. 

The quadratic map T can be rewritten as T = M^M\ by use of the product of two involutions 
Mi and M 2 , where M\ = Mf = 1@, § and, 

J I' = /-Ksin(#), 
Ml • \ 9' = -9, (mod 2tt) (2) 

M2 : \ 9' = -9 + 2vr/i - I 2 . (mod 2vr) (3) 
Moreover the mapping T has the symmetry S, i.e., ST = T^H, where 

S: \ 9' = 9 -tv. (mod 2vr) ^ 



2 The reconnection of twin-chains 



Since the twist function /^(/) is quadratic, the mapping T does not satisfy the twist condition 
(dfu/dl 7^ for V I), and the twin resonant chains appear on both sides of a shearless KAM 
curve g |, §. In Figs.l and 2, period -one and period-two twin-chains are shown for various 
values of K, with [i fixed. When we change the value of the parameters [x or K, the twin-chains 
exhibit topologically different structures; Fig. 1(a) for K = 0.05 shows that two twin-chains 
with period-one are generated on both sides of the shearless KAM curve, and Fig. 1(c) for 
K = 0.20 shows the saddle-node fusion of those two twin-chains. At the critical parameter 
value K = 0.1173906615, the separatrices of these twin-chains merge in the manner of the 
double homoclinicity(see Fig. 1(b)). For general cases of odd periodic twin-chains, the merging 
process occurs in the same way as described above. The parameter values at which the twin- 
chains merge are called reconnection thresholds. The reconnection thereshold for the period-one 
twin-chains, Krel, has been derived by Howard and Hohs by using the averaged Hamiltonian 
method as follows^, |4|: 

K^ c = \{2^fl\ (5) 

Although Eq.(5) was derived by the perturbation theory, the results explain the numerical data 
very wellQ. 

On the other hand, the twin-chains with even-periods reveal quite different features from twin- 
chains with odd-periods. Figure 2 is the case for the merging of the twin-chains with period-two. 
In the case of the even periodic twin-chains, the pair of the resonant chains appear with the 
same phase near the shearless KAM curve. The topological structure of the twin-chains is quite 
different from the odd periodic case. At a critical parameter value, the hyperbolic periodic 
points of the twin-chains merge(see Fig. 2(b)). We will call this state the reconnection for the 
twin-chains with even-periods. The reconnection threshold for the period-two 
been analytically derived as|4j 

^ = 2^2^ + 1) . (6) 

3 The symmetry of the quadratic map and the shearless KAM 
curves 

Hereafter we focus our attention on the KAM curves between twin-chains and derive the location 
of the shearless KAM curve. Denote the KAM curve which goes through the point (9,1) = (0, Iq) 
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as /C/ , i.e., 

£/o = {(*> !) = r " x o, x = (9 = 0, 1 = Jo), -oo < n < oo}. (7) 
The KAM curves are characterized by the rotation number R, 

1 N 

R = R(T, (0 O , /o)) = Um ^ " *n-i). ( 8 ) 

n=l 

which is a function of the mapping T and initial conditions (0o,Io). By fixing the value of 0$, 
the KAM curve is parametrized by the value of Jo- 
in the quadratic map, the KAM curve remains invariant under the involutions Mi, M 2 , which 
is shown as follows. By using the relation MfT n = T~ n Mj(j = 1, 2) for any integer n, 

MiT"x = T~ n Mix = T""x (9) 

and 

M 2 T n x = T-"M 2 x = T~ n (0 = 27T/x - I 2 , / = J ) = r^"" 1 )^, (10) 
where xo = (0,1) = (0,Iq)- Therefore we have 

MjKi = {(0,I)\(0,I) = MjT^xcXo = (0 = 0,1 = I ),-oo < n < 00} (11) 
= {(0,I)\(0,I) = T n x ,x = (0 = 0,1 = I ), -00 < n < cx)} (12) 
= /C /0 . (i = 1,2) 

The shearless KAM curve which passes through the point (0 = 0, 1 = Iq), say /C/*, is defined 
as the KAM curve whose rotation number R* satisfies the condition 

R(T, (0 = 0,1 = Jo)) + R* for V J / / *. (13) 

This definition is trivial in every part of the phase space if there exists only one shearless KAM 
curve in the phase space. 

Moreover, it is shown that the shearless KAM curve is invariant under the transformation S 
as follows. By using the relation ST m = T m S, we have 

R(T,S(0o,I o )) = R(T,(0 o ,I o )). (14) 

Therefore, if the shearless KAM curve goes through the point (9q, Io), it also goes through the 
point 5(^0,-^0)) i-e., the shearless KAM curve is an invariant set of S: 

SJC IS =1C IS . (15) 

We require that the shearless KAM curve )C[* crosses #-axis at two points, denoted by 
A : (0,1) = (0a,O) and B : (0,1) = (0b, 0) respectively. Since the shearless KAM curve is an 
invariant set of M2 and S, the following two equations hold: 

0B = -9a + 2^, (16) 

9 B = 9 A - it. (17) 

From the above equations, we have 

A = *{n + \), (18) 

9 b = k(v-\), (19) 

which are K independent. Equations (18) and (19) determine the exact location of the shearless 
KAM curves in the phase space. 
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4 The breakup of the shearless KAM curves 



In order to obtain the breakup threshold of shearless KAM curves, we numerically iterate the 
point (0,1) = (ir(fi + 1/2), 0) for each set of fi and K, and examine whether the motion is 
bounded or not. In cases when the motion is unbounded, the shearless KAM curve is con- 
sidered to have been destroyed. In our numerical calculations, the shearless KAM curve is 
considered to be broken up if the absolute value of I exceeds 2 during 10 5 iterations. 

The numerical results are shown in Fig. 3, where the motion is bounded in the black region, 
but the shearless KAM curve does not exist in the white region. This straightforward numerical 
method enables us to see the breakup condition of the shearless KAM curves. This method 
cannot detect the detailed threshold boundary for the breakup of the shearless KAM curves. 
However, the result will be much more refined by improving the breakup criterion in the numer- 
ical calculations. Actually, the breakup thresholds obtained in Fig. 3 are numerically confirmed 
by the determination of the rotation number for the nearby KAM curves. 

The breakup threshold of the shearless KAM curve characterized by the inverse golden mean 
I/7 rotation number has already been determined by del-Castillo-Negrete et al. using the 
Greene's residue method ||. In the case of the critical value for the destruction of the I/7 
shearless KAM curve, they obtained fj, c = —0.313951 and K c = 1.54156, which are indicated by 
the point (fj, c ,K c ) in Fig. 3. The shearless KAM curve for fj, c and K c is plotted in Fig. 4, which 
is obtained by numerically iterating the point [6 = 7r(^ c + 1/2),/ = 0). 

At the reconnection threshold of twin-chains, the rotation number of the shearless KAM 
curve, which exists between two twin-chains, is equal to the rotation number of the twin-chains. 
Thus, the reconnection threshold for the twin-chains with any period can be obtained by nu- 
merically calculating the rotation number of shearless KAM curves. In the cases of period-one 
and period-two, it is possible to derive the reconnection threshold analytically, as was shown 
previously. The P/Q lines plotted in Fig. 3 are numerically determined reconnection thresholds 
for twin-chains with the rational rotation numbers P/Q- In the quadratic map, the rotation 
number of a shearless KAM curve depends on the values of both fi and K. This is because K 
not only represents the strength of the non-integrable perturbation j7j, but also the change in 
the integrable part of the system. 

As shown in Fig. 3, the boundary of the breakup diagram reveals many sharp singularities on 
the lines of the reconnection thresholds with rational rotation numbers P/Q- Figure 5 shows a 
magnification of the breakup diagram around (fi c ,K c ), where fine singular structures are seen. 
We expect that much more fine structures, corresponding to the reconnection thresholds for 
twin-chains with higher periods, are immersed in the breakup diagram. 

5 Summary and discussion 

In this paper, we have determined the exact location of the shearless KAM curve in the quadratic 
map. Furthermore, we succeeded in obtaining the global structure of the breakup diagram of 
shearless KAM curves in the entire parameter space. The breakup diagram has many sharp 
singularities, and they are well explained by the reconnection thresholds for twin-chains. Fine 
structures as shown in Fig. 5 suggest that the breakup diagram might be a fractal. The theory 
which combines the breakup of shearless KAM curves with the reconnection of twin-chains 
remains to be established. 

As mentioned by del-Castillo-Negrete et al., a shearless KAM curve exhibits self-similar 
hierarchical structures for certain parameter values (H . Now that the exact location of a shearless 
KAM curve is obtained, it can be used to study the geometry of shearless KAM curves. The 
detailed analysis of the breakup diagram and the geometrical properties of shearless KAM 
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curves will be reported elsewhere. 
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Figure captions 

Fig.l. The merging process of the twin-chains with period-one for various values of K at fi = 
0.05. (a) if=0.05(below threshold). (b) i-C=0.1173906615(reconnection threshold), 
(c) i-C=0.20(above threshold). 

Fig.2. The merging process of the twin-chains with period-two for various values of K at 
H = -0.49. (a) K=0.4(below threshold), (b) K=0.50132565(reconnection threshold), 
(c) i-C=0.7(above threshold). 

Fig. 3. The breakup diagram of shearless KAM curves and the reconnection thresholds for 
the twin-chains with rational rotation numbers P/Q. The critical point (/x c = 
—0.313951, -fT c = 1.54156) obtained by del-Castillo-Negrete, Greene and Morrison is 
indicated by an arrow. 

Fig.4. The shearless KAM curve for (fi c ,K c ). 

Fig.5. A magnification of the breakup diagram around (fi c ,K c ). 
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Fig.l 




Fig.2 
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